Mesoscopic oscillator in U-shape with giant persistent current 
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A mesoscopic oscillator in U-shape has been proposed and studied. Making use of a magnetic flux 
together with a potential of confinement, the electron contained in the oscillator has been localized 
initially and an amount of energy has been thereby stored. Then a sudden cancellation of both the 
potential and the flux may cause an initial current which initiates a periodic motion of the electron 
from one end of the U-oscillator to the opposite end, and repeatedly. The period is adjustable. The 
current associated with the periodic motion can be tuned very strong (say, more than two orders 
larger than the current of the usual Aharonov-Bohm oscillation). Related theory and numerical 
results are presented. 



Due to the great progress in experiments, a few given 
number of electrons can be captured and confined in vari- 
ous ingenious artificial mesoscopic devices (say, quantum 
dots, wires, rings, and more complicated coupled dots, 
Mobius rings, etc.) [E 0, & i, 0, 0, 0, 0- These de- 
vices are basic elements in micro-industry. In developing 
micro-techniques a crucial point is the counting of time. 
In this paper a micro-oscillator in U-shape is proposed, 
which might work as a mesoscopic pendulum. 

It is recalled that, for classical motion, the crucial point 
of an oscillator is the storage of an amount of energy 
which can be transformed to kinetic energy later. For 
the classical pendulum, an amount of potential energy 
from gravity has been stored in advance, and will be 
transformed to the energy of swinging afterward. For a 
quantum oscillator, the crucial point is also how to store 
the energy. For this purpose, based on the idea suggested 
in the ref. [§], a device is proposed as sketched in FigQ] 
The device is a ring together with two arms (AB and 
DE), wherein a few free electrons are contained. Previ- 
ously, the arms are exactly blocked (say, by electrodes 
with high voltage), the ring is threaded by a magnetic 
flux <!>, and a strong potential V is applied so that the 
electron is localized in the bottom of the ring (close to 
point C in Fig. la). In this way, an amount of energy 
(which can be quite large) has been stored as shown later. 
Suddenly, both $ and V are cancelled, and the block on 
the arms is released. Instead, the upper half circle of 
the ring is blocked. With this sudden change, the ring 
is transformed to a U-oscillator (FigQ}}). The previously 
stored energy will motivate an oscillation of the electrons, 
as we shall see, from one end (A) to the other end (E). 
Dissipation is assumed to be negligible (namely, a super- 
conducting device). Then, the oscillation will proceed on 
in an exact periodic way. Related theories together with 
numerical results are as follows. 

For simplicity, it is assumed that only one free electron 
with an effective mass m* is contained in the device. Let 



FIG. 1: (Color online.) A sketch of the oscillator. The device 
in (a) is in fact a ring, where the two arms (AB and DE in 
dash line) are blocked. Whereas in (b) the upper half circle 
(in dash line) is blocked and the previous block on the arms 
is released. Thereby the ring is transformed to a U-pipe. 



the radius of the ring be R and the lengths of both arms 
be I. The ring and arms are very narrow so that they 
are quasi-one-dimensional. Previously, the Hamiltonian 
defined on the ring associated with FigQJi reads 



Hi, 
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V(6) 
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Where 9 is the azimuthal angle of the electron (the vector 
OO' pointing up in FigEQi has 6 = 0), G = j^z- $ 
is the magnetic flux in the unit <&o = hc/e. V{&) = if 
it — d < 6 < n + d, or = V Q otherwise. Where d measures 
the width of the confinement, V is a sufficiently large 
positive number so that the electron in the ground state 
tpg with energy E g > is strongly localized previously. 

In order to obtain ip g , Hi n n is diagonalized by using 
the set |fc) 



ik9 



/v2tt as basis functions, where k are 



integers ranging from 
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For numerical calculation, the range of k must be 
limited, and it was found that —40 < k < 40 is suf- 
ficient to provide accurate results (say, have at least 
four effective figures.). After the diagonalization E g and 
ip g = J2k Ck\k) can be obtained. 

Suddenly, the well is removed and V(9) becomes zero 
everywhere, the flux $ is also removed, and at the same 
time the block on the arms is released while the upper- 
half of the ring is exactly blocked. This leads to a change 
of the path of motion. Accordingly, the Hamiltonian is 
changed from Hi n u to 



H, 



evol 



ti 2 d 2 

2m* ds 2 
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defined in the U-pipe, where s is the distance of the elec- 
tron apart from A (say, when the electron locates at the 
lower-half ring, s = I + R(9 — tt/2)). After the change 
ij)g is no more an eigen-state of the new Hamiltonian, 
therefore the electron begins to evolve. 

The formal time-dependent solution of the new Hamil- 
tonian with the initial state ip g reads 



*(«,*) 



-iH evot t/h 
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This formal solution can be rewritten in an applicable 
form if we know all the eigen-states of H evo i . These eigen- 



states read simply \j) 



'2/Lsm(pjs), where L = 21 



Rtt is the total length of the U-pipe and s is ranged from 
zero to L, pj — jir/L and j is a positive nonzero integer. 
With them, Eq.(j4]) becomes 



*(M) = E< 



-iEjt/h 



(5) 



where Ej — h 2 p 2 /(2m*). The summation of j in Eq.((5|) 
is in principle from 1 to +oo. However, in numerical 
calculation j can be confined within a range (say, from 1 
to 80). 

From Eq.(f5]) the time-dependent density 
p(s,t) = |*(M)| 2 



unit of time. The time-dependent current defined from 
the conservation of mass reads 

J ( s > t) = i T,j,f cos(Jsn/L) sin(/s7r/L) j 



Im{e 



where only the imaginary part is contributed. In Eq.([7]) 
the unit of current is t~ 1 . From now on, we shall use r to 
measure the time. Obviously, both p and J are strictly 
periodic, the period of r is 2tt. 

Since ip g is localized in the bottom of the U-pipe, it is 
sufficient to carry out the integration involved in (j\ip g ) 
only in the domain tt/2 < < 3n/2. In this domain the 
eigenstates of H evo i can be rewritten as 



{-l)^-^' 2 ^2/L 
cos[jnR{6 - n)/L] 



(if j is odd) 



(8) 



li) 



/ 2/Lsm[jTrR(0-Tr)/L] (if j is even) (9) 



From Eqs.® and ([5]), obviously, |j) is symmetric with 
respect to 6 — ir if j is odd, or antisymmetric if j is even. 
Since the real (imaginary) part of |A;) is symmetric (anti- 
symmetric) with respect to ir, and since ip g is a superpo- 
sition of | fc) , the symmetries Eqs.® and ([9]) lead to a fact 
that (j\ip g ) would be a real number if j is odd, or an imag- 
inary number if j is even. Therefore, when j ±j' are even, 
the product {ipg\j'){j\ipg) is a real number, and the time- 
dependent factor e^^^i )0'+j ) T in Eq.© can be thereby 
rewritten as cos[(j - + j')r], where (j - + j') 
is an even integer I e . Alternatively, when j ± f is odd, 
the product (ipglj^ijl^Pg) is an imaginary number, and 
the time-dependent factor can be thereby rewritten as 
- sin[(j - + j')r] where (j - + j') is an odd 
integer I a . Since cos[/ e (7r/2 — t)] = cos[/ e (7r/2 + r)] and 
sin[/ (7r/2 — t)] = sin[/ (7r/2 + r)], we arrive at an im- 
portant feature of the evolution, namely, 



p(s,n/2-T) =p(s,tt/2 + t) 



(10) 



which implies a symmetry of time reflection with rc- 
Ejj' sin(js7r/L) sin(j's7r/L) (6) spe ct to tt/2. Furthermore, since sin\pj(L - s)] = 



where only the real part of the right hand side is con- 
tributed, r = t/t D , and t a = 2m* L 2 / {tt 2 Ti) is used as a 



(— 1)' J+1 svn(pjs) and cos[pj(L 
have 



-1)° cos(pjs), we 



p(L - S, 7T — t) = p(s, 7T + t) 



(11) 
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FIG. 2: (Color online.) The time-dependence of the aver- 
age position of the electron (s) and the average current J aV e 
plotted in a quarter of a period. The curves "1" to "4" have 
$ = 1.25, 3.25, 6.85, and 8.85, respectively. 

which implies a symmetry of time reflection with respect 
to 7r together with a spatial reflection with respect to 
L/2. 

Similarly, one can prove 

J(s,7r/2-r) = -J(s,it/2 + t) (12) 
J(L-s,ir-r) = J(s,tt + t) (13) 

Due to Eqs. (fT0|) and (fT2|) , the evolution in the duration 
[•7r/2,7r] is the time-reversal of that in [0, 7t/2]. Due to 
Eqs.fTTj) and (fl3|) . the evolution in [ir, 2ir] is the time- 
reversal of that in [0, it] together with a spatial reflection 
against L/2. Therefore, the whole evolution can be un- 
derstood if that in [0, tt/2] is clear. In what follows the 
study is restricted in [0, tt/2]. 

In order to have numerical results, let m* = 0.063me 
(for InGaAs), L = 2000nm, I = 400nm, R = 1200nm/7T, 
and V = 30G = 0.1243meV", d = tt/4. With these 
parameters, the time unit t = 4.408 x 1CP 10 sec. These 
choices are rather arbitrary. 

In order to have a general impression on the evolu- 
tion, we define the time-dependent average position of 
the electrons 

(a) = J P (s,t)s ds (14) 

and define the time-dependent average current as 

J a ve{r) = j- J J(s,t) ds (15) 

When <I> is given at four values, (s) and Jave(i~) against r 
are shown in FigO When r = 0, obviously (s) — L/2 (at 



FIG. 3: (Color online.) Spatial distribution of the densities 
p(s,r) and currents J(s,r) against s. The curves "1" to "4" 
have t = 0, vr/48, vr/24, and vr/12, respectively (these values 
are within the duration of oscillation DoS). The preset flux 
$ = 3.25 (the same for FigsO to |5J). The unit of current is 
t^ 1 , where t a = 2m* L 2 /(iy 2 Tl) is the unit of time used in this 
paper. 



the bottom) as shown in Fig[2^. However, the current is 
not zero initially (J a ue(0) 7^ 0). This is shown in FigJ^b, 
where a larger $ leads to a more negative J a ue(0). 

Let us define E evo i = {ijj g \H evo i \ip g ) which is the energy 
contained in ip g after the transformation of the Hamilto- 
nian. It was found that, when <I> increases from 1.25 to 
8.85, E evo i increases from 0.0152 to 0.333meV. Obvi- 
ously, a larger E evo i will lead to a stronger initial current 
which motivates the evolution afterward. 

In FigHthe duration [0, tt/2] can be roughly divided 
into two, [0, 7t/6] and [77/6, 7r/2]. In the former (s) goes 
toward the two ends alternately. Accordingly, J ave ap- 
pears as negative and positive alternately. Obviously, 
this implies that the electron oscillates end-to-end repeat- 
edly. When $ is larger, the first minimum of (s) will shift 
down and left as shown in Fig|2^. It implies that, once 
the evolution begins, the electron will be closer to A in 
a shorter time. Furthermore, a larger $ will cause more 
rounds of oscillation taking place in the duration [0, 7t/6]. 

In the duration [tt/Q, tt/2), (s) remains to be close to 
L/2, and | J a ve\ becomes small. It implies that the proba- 
bilities of the electron staying at the left and right sides of 
the bottom are nearly equal, and both negative and pos- 
itive currents might appear in the path simultaneously 
(this leads to a cancellation and thereby a smaller | J a ve\)- 
The above two durations, for simplicity, are called dura- 
tion of oscillation (DoS) and duration of cruise (DoC). 

To study in detail the evolution in the DoS, the spa- 



4 




0.0 2 0.4 0.6 0.8 1 



s/L 

FIG. 4: (Color online.) The same as Fig[3]but in the duration 
of cruise DoC [tt/6, n/2]. The curves "1" to "3" have r = 
0.2857T, 0.3547T, and 0.4587T, respectively. 

tial distribution of p(s,r) is plotted in FigJ31 where r is 
given at four values in the early stage of evolution. When 
t = 0, the curve " 1" of FigJSJi shows the initial localiza- 
tion, and " 1" of FiglSJ} shows the strong negative initial 
current created via the sudden change of the Hamilto- 
nian. Due to the negative current, the peak of p shifts 
left rapidly as shown by "2" of Fig[3^. Correspondingly, 
the distribution of J shifts also left as shown by "2" of 
FiglSb- Afterward, the peak of p keeps going left and will 
be close to A as shown by "3" of Fig[3K However, "3" 
of Fig|3]3 is mainly positive. It implies that the direction 
of motion has already been reversed and the electron be- 
gins going toward the other end E. A little time later, the 
density is partially close to E as shown by " 4" of FigJHK- 
Meanwhile, the current is positive throughout the path 
as shown by "4" of FigJ3jD. Accordingly, the distribution 
of p as a whole is going right. FigJ3] shows only the first 
round of oscillation in the DoS when $ = 3.25. 

The maximal initial current J(L/2,0) (e.g., the dip 
of "1" of Figl3j3) was found to be nearly linearly pro- 
portional to <!> (with our parameters, J(L/2,0) w 
— 5.69$/t ). Incidentally, if $ < 0, the initial current 
would be positive and therefore the direction of motion 
would be reversed. 

In the DoC both p and J, in general, are widely dis- 
tributed along the path from A to E with numbers of 
peaks and dips as shown in Figj4] In this duration the 
current may be positive (going right) somewhere and neg- 
ative (going left) elsewhere. If J = and J > take 
place at s — s a , there would be a source at s from where 
the current flows out to both sides. Whereas if J = and 
■S-J < at s , there would be a leak to where the cur- 
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FIG. 5: (Color online.) The same as Fig|3]but the curves "1" 
to "3" have r = 7r/4 — 7r/96, 7t/4, an d 7r/4+7r/96, respectively. 

rents flow in from both sides. Both sources and leaks are 
found in the path. The classical picture of motion of the 
electron in the DoC is not clear, it seems to be chaotic. 

However, in the DoC, there is a noticeable instant. 
When t — tt/A, we have (s) = L/2 and J ave — disre- 
garding how $ is as shown in Fig[2j At this instant p is 
no more widely distributed but concentrated close to the 
two ends and the bottom C as shown by "2" of FigfS^ 
(where the peak at the bottom is much higher). When r 
is a little earlier and later than tt/4, p and J are shown, 
respectively, by "1" and "3" of Figdji and [5b. They to- 
gether demonstrate how the density is concentrated into 
three peaks and spread out afterward. The abscissa of 
FigJSJa has been divided into four regions. The curve 
" 1" of FigfSb is negative in region I. Thus the density is 
pushed left resulting in forming the peak of p close to A. 
Besides, " 1" is positive in II but negative in III. Thus p 
is pushed from both sides of C to the bottom resulting 
in forming the highest peak. Furthermore, the positive 
current in IV leads to the peak close to E. Exactly at 
the instant r — ir/4, the current is zero as shown by "2" 
of FigJSJa. It implies that the system is static instanta- 
neously. Afterward, the current increases rapidly but in 
reverse direction as shown by "3" of FigfSjD. This leads 
to the disappearance of the three peaks as shown by " 3" 
of Figj5^. The existence of instants wherein the system is 
instantaneously static is not at all surprising, this might 
happen in various types of oscillation (say, in classical 
pendulum) . 

When r is close to 7r/2, the density is once again con- 
centrated in the bottom as shown by " 2" of Fig[6K, where 
the main peak is very sharp and high. At the same time, 
the associated current (" 2" of FigEJa) is exactly zero (this 
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FIG. 6: (Color online.) The same as FigObut the curves "1" 
to "3" have r = 7r/2 — 7r/96, 7r/2, and 7r/2 + 7r/96, respectively. 

is obvious from Eq . (fl~2"|) ) . Thus the system becomes once 
again static instantaneously. "1" and "3" of 6a describe 
the densities a little earlier and later, respectively, than 
7r/2. They overlap exactly due to the symmetry given in 
Eq.Cni. "1" and "3" of Figib are exactly opposite to 
each other due to Eq. (fT2")) . In fact, they cause a rapid 
gathering and, successively, a rapid extension of p. In 
particular, the exact symmetry appearing in Figj6j} im- 
plies that the evolution has arrives at a turning point. 
Afterward, the evolution will proceed exactly reversely as 
demonstrated by Eqs. UU\) and Alty) . For examples, the 
density will be re-gathered into three peaks exactly as 
"2" of FigGJi when r = 37r/4, and will re- visit A exactly 
as "3" of 3a when r = tt — 7r/24, etc. When r = tt, p is 
exactly as " 1" of 3a, while J is exactly opposite to " 1" of 
FigJHb. Thus, instead of going left, the peak of p at the 
bottom goes right at r = tt due to the positive J. 
It is noted that from Eqs. (fT0]) and (fTTj) . we have 

p(L - s, r) = p(s, tt + t) (16) 

From (|T2"]) and (fTB"]) , we have 

J(L-s,t)=-J(s,tt + t) (17) 

In the duration [tt, 2tt], Eqs. (fit)]) and (JTTJ) together imply 
that the evolution is just a repeat of that from to tt, but 
with an interchange s < — > L — s (namely, an interchange 
of the A and E ends). Thus, the evolution in the whole 
period 2tt is completely clear, and it goes on again and 
again periodically. 

Now let us compare the current in the U-oscillator 
with the famous Aharonov-Bohm (A-B) persistent cur- 
rent of an electron on a ring. The maximal A-B cur- 
rent of the ground state is Ja-b = h / '(Anm* ' R 2 ) . Thus, 
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FIG. 7: (Color online.) p(s,r) and J(s,r) with $ = 8. The 
curves "1" to "4" have r = 0, vr/96, tt/48, and tt/32, respec- 
tively. 

i^ 1 = 2tt(L — 2l) 2 /L 2 J A-B- For our parameters, we 
found that J(i/2,Q) = -12.87$ J A -b- Thus, giant cur- 
rent can be obtained if $ is sufficiently strong (say, if 
<f> w 10, J(L/2,0) is two orders stronger than the A-B 
current). 

Since t is proportional to m* L 2 , a longer path and/or 
a heavier effective mass will lead to a longer period. The 
preset $ as a motivity is crucial to the oscillator. Its effect 
has been shown in Fig [2] and is further shown in Fig [7] 
with $ = 8 (to be compared with Fig.3 with $ = 3.25). 
Comparing the distributions of "3" of Fig[7H and FigEK, 
the former is closer to A and occurs much earlier (the 
electron rushes to A more rapidly). While $ is a very 
sensitive factor, the parameters of the preset potential is 
less sensitive to the evolution. 

In summary, a U-oscillator has been proposed and 
studied. The motion of electron in the device is strictly 
periodic and adjustable. The following points are re- 
minded 

(i) When t a is used as the unit of time, the period 
is 2tt. The evolution in the second quarter of a period 
is the time-reversal of that of the first quarter, and the 
evolution in the second half is the same as the first half 
but with the interchange s < — > L — s, namely, a spatial 
reflection against L/2. 

(ii) Each quarter can be divided into two durations, 
namely, DoS and DoC. The classical picture in the former 
is clear (a few rounds of oscillation appear), but not clear 
in the latter. When $ is larger, more rounds of oscillation 
will appear in the DoS. 

(iii) Let the time associated with the lowest minimum 
of Fig. 2a be T en d- Then, in the first half period, the elec- 
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FIG. 8: The time-dependence of (s) in a whole period. $ = 
3.25 (a) and 8 (b). 

tron is very close to A twice at T enc i and tt — T enc i, respec- 
tively. Whereas in the second half period, the electron is 
very close to E also twice at n+T enc i and 2tt— T enc i, respec- 
tively. Such a close contact of the electron with the two 
ends in a whole period is shown in FigjEJ where the ap- 
pearance of rounds of oscillation in the DoS is also clear. 
On the contrary with classical pendulum, it is clear from 
FigJS] that is explicitly the end-to-end oscillation occurs 
only if r is close to In (J is a positive integer). 

(iv) Each time when t — In/A , the system is instan- 
taneously static, namely, J(s,Itr/4) = for all s. Mean- 
while, the density would be mostly concentrated in the 
neighborhood of C, namely, the bottom of the U-pipe. 



(v) The introduction and the sudden removal of $ is 
a crucial point. This leads to the sudden creation of the 
initial current, which motivates the evolution afterward. 
In particular, when <E> is sufficiently large, giant periodic 
current (two or more orders stronger than the A-B cur- 
rent) can be obtained. 

The device might work as a quantum pendulum. It 
can be generalized to include a group of localized elec- 
trons initially. This case deserves to be further studied. 
The idea proposed in this paper might be useful in micro 
technology. 
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